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Abstract: We study the diffractive electroproduction of quarkonia from 
quark- ant iquark-gluon states in the photon wave function. We show that these states 
contribute to the leading-power and leading-logarithm level. We suggest the mea- 
surement of J / ip production via inelastic diffraction to study color-transparency and 
color-opacity effects in the diffractive gluon distribution. 
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1. Introduction 

J/ijj mesons are produced copiously via hard-diffractive reactions in lepton-hadron 
collisions at high energies. In these processes, the initial hadron emerges intact in 
the final state and the meson is produced via a short-distance interaction, involving 
momentum transfers of the order of the heavy quark mass or, for highly off-shell 
photons, the photon virtuality. 

The elastic process 7 + A — > V + A', where A and A' are the initial and final 
hadrons and V is the quarkonium, has been studied extensively at HERA [1, 2, 
3]. This process is computable in QCD [4, 5, 6] in the sense that a short-distance 
coefficient, calculable as a power expansion in a s , can be factorized [7] from a parton 
density describing the long-distance dynamics of the hadronic state A and a wave 
function describing the long-distance dynamics of the produced bound state. A 
typical leading-order graph is shown in Fig. 1. The parton density probed by this 
process is a nondiagonal parton density [7, 8, 9, 10]: 

4 a (Z,x F ,t,ri = ^+ / dye^iy- (A'\G a (0) +j G a (0,y-,0) +j \A) , (1.1) 

where p\ is the lightcone plus momentum of hadron A, x P = 1 — p\ijp\i t = 
{Pa — Pa') 2 - The operator G(0,y~ , 0) is the gluon field strength multiplied by a 
path-ordered exponential of the line integral of the color potential from y~ to 00. 
The ultraviolet divergences from the operator product are renormalized at the scale /x. 

In this paper we are concerned with the production of quarkonia by inelastic 
diffraction, 7 + A — > V + X + A', a process that has not been considered so far 
but could be measured at HERA and at future lepton-proton and lepton-nucleus 
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colliders. This process probes a different parton density than (1.1): see Eq. (2.1) 
ahead. As we will see, for sufficiently large values of the hard-scattering scale its 
cross section is enhanced by two powers of the hard scale compared to the case of 
Fig. 1. 

In the rest frame of the target (see Sec. 2), the elastic production is dominated by 
fluctuations of the incoming photon into quark-antiquark states, while the process we 
consider in this paper is dominated by fluctuations into quark-antiquark-gluon states. 
The main interest in measuring this process is to observe effects of the coupling of 
these states to the target's soft color field. The color dipole interaction is stronger 
for QQg states than for QQ states. We discuss consequences of this in Sec. 5. 

The contents of the paper is as follows. We present the physical picture of the 
process in Sec. 2. This discussion uses the target rest frame and gives an "s-channel 
view" of the process, based on the evolution of the system into which the incoming 
photon dissociates. Then we relate the evolution of this system to diffractive matrix 
elements of gluon field operators, and in Sees. 3 and 4 we give factorization formu- 
las for the quarkonium cross sections. This formulation corresponds to a "t-channel 
view" of the process. We use these results to contrast the scaling behaviors of the 
elastic and inelastic contributions. Sec. 5 contains a concluding discussion, focus- 
ing on the possibility that diffractive interactions of QQg states with the target be 
dominated by small transparency lengths, and on the implications of this for the 
diffractive gluon distribution. The analysis of this paper applies to both J ftp and T 
mesons. Numerical evaluations of the cross sections are left to future work. 

2. J/V 7 production by QQg states: the s-channel point of view 

In this section we work in the rest frame of the target and present the physical 
picture for the production of quarkonia from quark-antiquark-gluon states in the 
photon wave function. 

As is well known, the spacetime structure [11, 12] of hard-diffractive processes in 
electroproduction looks especially simple in the target rest frame. In this frame, at 
lightcone times far in the past a highly energetic photon dissociates into a partonic 
system. Much later, this system interacts with the color field of the target, mainly 
through its color dipole moment [12]. The physics of the parton densities is re- 
expressed in this frame in terms of the dipole-target scattering (see reviews in [13, 
14]). The elastic production of quarkonia (Fig. 1) has been analyzed in this context in 
Refs. [4, 5, 6]. This process is dominated by the contribution of the quark-antiquark 
component in the photon wave function. 

Fig. 2 depicts a typical inelastic contribution to the diffractive production of 
quarkonia. In the target rest frame, a quark-antiquark-gluon state is created by the 
photon well outside of the target hadron. This system has large minus momentum 
and travels a long distance before interacting with the target's color field. (This 
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interaction is schematically represented in Fig. 2 as two-gluon exchange. In general, 
it includes any number of soft gluons and is nonperturbative. Note also that con- 
tributions in which k\ and k 2 attach to the quark lines in the upper subgraph are 
to be considered together with the contribution depicted in Fig. 2. One needs to 
include them to construct gauge-invariant quantities.) Over long times, the quark 
and antiquark bind into a quarkonium. The emission of two gluons into the final 
state in Fig. 2 is required by color conservation. The interaction of the fast-moving 
parton system with the target's field leaves hadron A almost intact, and produces a 
final state consisting of the quarkonium plus jets plus the diffracted hadron. 




Figure 1: A typical graph for the elastic process 7 + A — > V + A'. 




Figure 2: QQg contribution to diffractive quarkonium production. 

In the long interval between the formation of the parton system and the inter- 
action with the target's field, gluons can be radiated, possibly at small lightcone 
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separations. These graphs contribute to evolution and radiative corrections. We will 
see later the result of including them. 

Consider now the upper subgraph in Fig. 2, in which k\ and k 2 couple to the 
QQg system. This subgraph is dominated by the region in which the transverse 
momentum |p| is of order M (with M the heavy meson mass) and much larger than 
the gluon transverse momenta |k| and |k'|. This region contributes a hard-scattering 
squared amplitude of order aa 2 s1 that is, of the same order as the hard scattering in 
Fig. 1. By power counting in the coupling, the contributions in Figs. 1 and 2 are 
of the same size. Later on we will see that the contribution in Fig. 2 dominates the 
contribution in Fig. 1 by power counting in the hard-scattering scale. 

The observation above indicates that in the dominant region the gluon k goes 
near the mass shell. The process becomes sensitive to long distances, and involves 
the nonperturbative distribution of gluons in the hadron A, subject to the condition 
that the hadron is diffractively scattered. This distribution is depicted in Fig. 3. It 
is defined in terms of the gluon operator G introduced in Eq. (1.1) by the matrix 
elements 



df, 



dx F dt ' ' ' (4ir) 3 f3x F pj[ 
x(A\G a (0) +3 \A',X)(A',X\G a (0,y-,O)^\A) , (2.1) 

with f3x F the fraction of the hadron's plus momentum carried by the gluon. These 
matrix elements have been studied in detail in [15, 16]. The vertex at the top of 
Fig. 3 represents the gluon operator in Eq. (2.1), including a color-octet eikonal 
line defined by the path-ordered exponential in G [15, 16]. The eikonal line can be 
interpreted physically as a stand-in for the quark-antiquark pair that was produced 
by the electromagnetic current but is too compact to be resolved by the hadron's 
color field. Attachments of gluons to this line are incorporated as explained in [15]. 





Figure 3: Diffractive matrix element for the gluon operator (Eq. (2.1)). The graph on the 
right hand side corresponds to the contribution depicted in Fig. 2. 
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3. The t-channel point of view 



The precise relation of the vector meson cross section with the matrix elements (2.1) 
is given by a factorization formula. This can be obtained by applying to the graph 
of Fig. 2 standard methods [17, 18] for determining the leading-power asymptotics 
in the limit of a large hard-scattering scale. 

The leading integration regions in momentum space involve (see Fig. 4) i) a beam 
jet in the direction of hadron A, including the final hadron A'; ii) a hard subgraph, 
in which all virtual lines are off shell by order M; iii) final states jets, including the 
vector meson, that are not in the direction of A; iv) a soft subgraph, consisting of 
soft gluons that link the beam jet with the hard subgraph and final state jets. 

The method [17, 18] instructs us to a) use Ward identities to factor out soft gluon 
subgraphs from final-state and beam jets, and b) sum over real and virtual graphs 
to show that all soft gluon exchanges (dashed lines in Fig. 4) cancel in the cross 
section. This can be accomplished in a standard manner in the problem at hand, 
except for the treatment of soft gluons coupling to the beam jet in the plus direction: 
owing to the diffractive requirement on the final state, there are both initial-state 
and final-state couplings, giving poles on both sides of the integration contour in the 
complex plane of the gluon minus momentum r - . This prevents us from deforming 
the r~ contour to the region where the soft approximation is valid. 




Figure 4: Factorization of the inelastic diffraction process. Soft gluon exchanges (dashed 
lines) cancel. 

This obstruction, however, can be overcome in the same way as in diffractive 
deep inelastic scattering [19], by going to the complex plane of the plus momentum 
component r + , and showing that in this plane the contour is not trapped by singu- 
larities in the beam jet subgraph. Then the cancellation of soft gluon exchanges goes 
through, and factorization follows. To keep the notation simpler we first consider 
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the case Q 2 = 0, limiting ourselves to the direct-photon part of photoproduction. 
(We comment on the resolved-photon part below.) The factorization formula for the 
•y + A^V + X + A' cross section reads 

^-^(M 2 ,W 2 ,x P ,t) = £ f df3-%(f3,x P ,t,^)a a (-^,M 2 ,^) 

(3-1) 

with W 2 = (pa +P-y) 2 , w = M 2 /W 2 . The coefficient a a describes the meson produc- 
tion subprocess represented by the upper subgraph in Fig. 4. In the nonrelativistic 
expansion for the bound state, a a factorizes into the product of a perturbative factor, 
describing the production of the quark- ant iquark pair, and a nonperturbative factor, 
describing the binding of the quark and antiquark into a quarkonium. In this paper 
we consider only the lowest order in this expansion, corresponding to color-singlet 
quark- ant iquark states [20]. We comment below on the possible role of color-octet 
states. 

Both a a and f a in Eq. (3.1) depend on the factorization scale ji. As noted 
earlier, the evolution in // arises, via renormalization of the ultraviolet divergences, 
from gluons being radiated from the top subgraph in Fig. 2. We refer the reader to 
Refs. [15, 16, 19, 7] for detailed treatments of evolution in hard-diffractive processes. 

To the leading order in a s , only the term with a = g enters in the sum in 
Eq. (3.1). This is the term depicted in Figs. 2,4. Compare this with the case of the 
diffractive structure function F 2 , where a — g enters to the next-to-leading order, and 
of diffractive jet production, where both a = g and a = q enter to the leading order. 
(See, e.g., [21, 22] for early treatments of gluon radiation in these processes.) This 
makes the quarkonium cross section an especially important experimental probe of 
the physics of large gluon densities in diffraction. 

The coefficient a g that relates the quarkonium cross section to the diffractive 
gluon density can be obtained from the perturbative calculations of [20] and reads 



32 a e 2 Q a s V) |^(0)| 2 



3tt M 2 M 3 

i r i 1 1 1 

(3.2) 



1 +LU 1-UJ ,1 

1 + M V? + 2UJ 111 - 

l-UJ (l+U)) 2 UJ 



(1+CcO 3 (1-C^ 2 



Here ip(0) is the quarkonium wave function at the origin. Eq. (3.2) vanishes like 
uj for large energies {uo — > 0) and like (1 — uo) near threshold (uo — > 1). In these 
boundary regions higher order corrections [23] to the coefficient become particularly 
important. Notice the behavior of the cross section (3.1) in the hard scale M 2 . Apart 
from the dimensionless factor \ip(0)\ 2 /M 3 , coming from the nonrelativistic treatment 
of the bound state, da/dt goes like 1/M 2 . This is to be contrasted with the behavior 
1/M 4 [4, 6] of the elastic process (Fig. 1). That is, for sufficiently large M the cross 
section is significantly larger than its elastic part. 



6 



As noted above, this treatment does not include terms from color-octet quark- 
antiquark states. Color-octet terms can in general be large, because, although sup- 
pressed by powers of the heavy quark velocity, they are enhanced by a power of 
l/a s [24]. In this respect, they are more important in the inelastic-diffraction pro- 
cess than in the elastic case. Observe however that if the final state is diffractive the 
contribution of color-octet states does not necessarily follow the simple factorized 
form that is commonly used in applications of the nonrelativistic expansion. The 
soft gluons emitted in the transition of the color-octet quark pair to quarkonium 
may couple to the diffracted final state. These soft gluons are likely to reduce the 
probability for producing a quarkonium via a color-octet pair while leaving hadron 
A intact. The net effect is that color-octet contributions should be suppressed in 
diffractive production compared to expectations valid in the fully inclusive case. 

The resolved-photon part of photoproduction cannot be described by a factor- 
ization formula of the kind (3.1), because the possibility of color being exchanged 
with the photon beam jet spoils the cancellation of the soft gluon subgraphs. Since 
the direct and resolved parts cannot be separated in general beyond the leading 
logarithms, this limits the applicability of factorization in the Q 2 = case. 

Note however that, if we estimated the sign and order of magnitude of the soft 
contributions from diffractive hadron-hadron data (see, e.g., [25]), this would suggest 
that the relative importance of the resolved part compared to the direct part is 
smaller in diffractive photoproduction than in inclusive photoproduction, perhaps 
by a factor of several. This, combined with the use of kinematic cuts designed to 
suppress the resolved process, may allow one to make a meaningful test of the direct 
process. 



4. Q 2 ^ case 

The leptoproduction case Q 2 ^ is conceptually simpler than photoproduction, 
because it does not involve the processes with soft color exchange discussed at the 
end of the previous section. In the leptoproduction case a result of the kind (3.1) 
holds. More precisely, both the longitudinal and the transverse cross sections satisfy 
a factorization formula in terms of the diffractive matrix elements (2.1): 

<laU "" l " ] (M 2 ,S,Q 2 ,y,x P ,t) = ^ (4-1) 



dx P dt dQ 2 dy " w yQ 

x 



l-y+ y 4)^r(M 2 ,Q 2 ,yS,x P ,t) - t--^-{M 2 ,Q 2 ,yS,x P ,t) 



^ =T,Jd(3^j t ((3,x P ,t,f, 2 )a ha (^,M 2 ,Q 2 ,f, 2 ) (i = 2,L). (4.2) 



dx P dt a ' dx P dt ' (3x, 
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Here, in standard notation, o-i is the sum of the transverse and longitudinal cross 
sections, a/5* is the lepton-hadron center-of-mass energy, y is the lepton energy loss 
in the target rest frame, and uj' = (M 2 + Q 2 )/{yS). 

Note that in the case of the elastic process (Fig. 1), in contrast, a factorization 
formula exists so far only for the longitudinal polarization [7]. 

The coefficients a ija generalize Eq. (3.2) to Q 2 ^ 0. Consider the power counting 
at large Q. (By comparison, recall [5, 6] that the elastic process of Fig. 1 scales like 
1/Q 6 at large Q.) For Q > M, a factor of 1/Q 2 replaces the 1/M 2 of Eq. (3.2). 
Fig. 2 has a suppression factor M 2 /Q 2 for producing a quarkonium of mass M when 
the momentum flowing in the quark-antiquark subgraph is of order Q. However, note 
that for Q ^> M Fig. 2 is no longer dominant, as contributions from the intrinsic 
heavy-quark distribution and from parton fragmentation into quarkonium scale like 
1/Q 2 . We therefore suggest considering the region of Q of order M, or smaller than 
M. In this region, quark-antiquark-gluon states give a contribution to diffractive 
quarkonium production that is both leading in powers of the hard scale and leading 
in the running coupling. 

Detailed studies are warranted to assess the numerical importance of this con- 
tribution at HERA and future eA colliders. 

5. Discussion 

We have seen that the electroproduction of quarkonia by inelastic diffraction mea- 
sures, through the factorization formulas of Sees. 3 and 4, the diffractive gluon distri- 
bution (2.1). For Q of order M, or smaller than M, this process provides a remark- 
able probe of quark-antiquark-gluon states in the photon lightcone wave function. 
Recall that, in contrast, the production of quarkonia by elastic diffraction couples 
predominantly to quark-antiquark states. 

The differences between the dynamics of QQ and QQg states scattering off the 
hadron's soft color field have been investigated extensively in the color dipole approx- 
imation (see, e.g., [12, 13, 14], and references therein). Color-octet dipoles interact 
more strongly with the target. As a result, the probability for the parton system 
to get through the hadron without breaking it up dies off faster as the dipole size 
increases, and the hadron starts to look opaque earlier. 

Based on the analysis of these effects, a physical picture for the diffractive par- 
ton distributions was suggested in [26] . One of the implications of the onset of color 
opacity is that, if the maximum parton separation 1/n for which the hadron is trans- 
parent is sufficiently small, a perturbative calculation [27] of the /3 dependence of the 
diffractive distributions should apply. Such a calculation is incorporated in [26] and 
in different models [28, 29] for the diffractive distributions. The overall normaliza- 
tion of the distributions is inversely proportional to the square of the transparency 
length 1/k (in units of the hadron radius r^), while their shape in (5 is given by 
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perturbatively computable functions ip: 

f a oc (N 2 - 1) K y A <p g (P) , f q ^N c K 2 q r\ <p q (P) . (5.1) 

In [26] predictions based on these distributions were compared with the HERA 
data [30] for the diffractive structure function F 2 . The comparison indicates that 
transparency lengths are smaller for QQg states than for QQ states. The F 2 data 
are consistent with K g /K q ~ Ca/Cf. Note that, together with the color factors 
explicitly shown in (5.1), this value gives back the result [27] for the ratio of the 
diffractive gluon and quark distributions in a small dipole, 

f 9 JC C\{Nl - 1) _ 27 
/, C 2 F N C 2 

thus supporting the picture based on the dominance of small parton separations. 
But a probe that, unlike F 2 , couples directly to gluons will give a much better way 
of exploring this picture. The vector meson cross sections considered in this paper, 
Eqs. (3.1) and (4.2), will allow one to study precisely the relation between the large 
size of the diffractive gluon distribution in a proton and the dominance of small 
transparency lengths. 

Note that this study may also be critical for deep inelastic scattering on nuclei, if 
information extracted from diffractive processes is to be used to compare approaches 
to nuclear shadowing [31]. 

It is worth commenting on the x P dependence. In conventional fits to diffractive 
data, a simple factorizing form of the x P and f3 dependence is often assumed. From 
the QCD point of view, this is an ansatz, not a theory result. In particular, notice 
that the diffractive gluon and quark distributions at low /x need not have the same x P 
dependence. Therefore, even if a simple factorizing form x p ~ 2a is taken, it would be 
interesting to analyze the data using two different behaviors for the gluon and quark 
distributions, parameterized by a g and a q . The heavy meson cross section discussed 
in this paper would provide a good way of determining a g . 
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